Evolution of entanglement for quantum mixed states 
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A simple relation is introduced for concurrence to describe how much the entanglement of bipartite 
system is at least left if either (or both) subsystem undergoes an arbitrary physical process. This 
provides a lower bound for concurrence of mixed states (pure states are included) in contrast to the 
upper bound given by Konrad et al [Nature Physics 4, 99 (2008)]. Our results are also suitable for 
a general high dimensional bipartite quantum systems. 
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I. INTRODUCTION 

Quantum entanglement is an important physical re- 
source in quantum information and computation tasks 
such as quantum teleportation [1], quantum key distri- 
bution [2], quantum computation [3] and so on. In real- 
istic quantum information processing, entanglement has 
to be prepared or distributed beforehand by two distant 
parties in which one or more physical systems have to be 
transmitted by a quantum channel. However, unlike clas- 
sical systems, quantum systems are usually fragile. It is 
inevitable that environment (channel) will influence the 
systems of interest more or less and induce decoherence 
because of the interaction with the systems, so that en- 
tanglement is destroyed to some extent before use [4-6]. 
It is a key task to evaluate the shared entanglement after 
the influence of environment. 

In usual one has to deduce the time evolution of entan- 
glement of the composite system from the time evolution 
of the quantum state under consideration [7], when the 
subsystems of the composite quantum system undergo a 
physical process. That is to say, for the different poten- 
tial initial states one has to repeat the same procedure 
every time. Quite recently, an important step has been 
taken by Konrad et al [8] who provided an explicit evolu- 
tion equation of quantum entanglement quantified by the 
remarkable concurrence [9] for bipartite quantum state 
of qubits. It has been shown that given any one-sided 
quantum channel, the concurrence of output state corre- 
sponding to any initial pure input state of interests can 
always be cquivalently obtained by the product of the 
concurrence of input state and that of the output state 
with the maximally entangled state as an input state. 
However, for two-sided quantum channel or the initial 
mixed states, the product of the two concurrences only 
provides an upper bound for the concurrence of inter- 
ests. It is obviously important to find a lower bound of 
the concurrence in order to well grasp the entanglement 
of output states. 

In fact, quantum mechanics and quantum information 
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processing (QIP) are not constrained to pure states as 
well as 2 (g) 2 dimensional quantum systems. Because of 
imperfect experiments and inevitable disturbance of envi- 
ronments, mixed (entangled) states are ubiquitous, which 
is inevitable for QIP to cope with. Furthermore, on the 
one hand, when we face the quantum features, especially 
entanglement, of meso- even macroscopic quantum sys- 
tems, such as Bose-Einstein condensates [10]. we have 
to deal with high dimensional density matrices. On the 
other hand, it has been shown that some QIP tasks based 
on high dimensional entangled states are more efficient 
than those of qubits. For example, cryptographic proto- 
cols are more secure based on quantum channel of qutrits 
[11-13]. Teleportation can be implemented in faith even 
though a non-maximally entangled quantum channel is 
shared [14]. Therefore, for these large systems, it is neces- 
sary to investigate the entanglement of high dimensional 
quantum systems (multi-parties included). 

However, quantification of entanglement, as a precon- 
dition of studying entanglement, is generally a hard prob- 
lem which does not only lie in the poor practicability for 
high dimensional quantum systems [15,16] but also the 
nonlinearity on density matrices for usual entanglement 
measures [17,18]. Therefore, it is often suggested to de- 
rive a lower bound to evaluate the entanglement (Entan- 
glement can be better evaluated, if both upper and lower 
bounds are given). In this paper, we consider the evo- 
lution of entanglement with either (or both) subsystem 
undergoing an arbitrary quantum channel. With concur- 
rence as entanglement measure, we find a lower bound of 
the concurrence of any output states for a given quantum 
channel based on the evolution of a probe state as input 
states. It is shown that our lower bound is not restricted 
to the bipartite systems of qubits. In particular, for one- 
sided quantum channel, our lower bound has a concise 
form. Furthermore, we also show that it is not necessary 
to choose the maximally entangled state as the probe 
states. The paper is organized as follows. In Sec. II, we 
give a lower bound for the concurrence of bipartite quan- 
tum systems; In Sec. Ill, we show that the lower bound 
can be obtained in terms of the initial input state and 
the output state with the probe state as initial states; 
The conclusion is drawn finally. 
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II. THE LOWER BOUND FOR CONCURRENCE as 



An (Aq ® N 2 )- dimensional bipartite quantum pure 
state can be written as 



Afi-l N 2 -l 

\^)ab = E E ^« ly) 

j=0 j=0 



(1) 



where denotes the computational basis and Aq x N 2 
matrix ip (i.e. IV^ab without \) AB ) represents the 
matrix notation [19,20] of \^) ab with matrix element 



1. Consider the 



iVi-l N 2 -l 
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Schmidt decomposition, IV^ab can a l so he given by 
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(2) 
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with A, being real in decreasing order and R = 
min{Aq, N 2 }. The concurrence of \4>) AB can he defined 
[21,22] as 



CMab) = 
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The equivalence of eq. (3) and eq. (4) lies in that quan- 
tum pure states are related to its Schmidt decomposition 
by local unitary transformations which do not contribute 
to concurrence. In particular, if N\ — N 2 = 2, eq. (4) 
can be reduced to C(\ip) AB ) = 2AiA 2 = 2 |det (-0) | , which 
will be used later. With the definitions given by eqs. (3) 
and (4), we obtain the following theorem. 

Theorem 1. For any (Aq <g> N 2 )- dimensional bipar- 
tite quantum pure state 



(5) 



with £ denoting the set of (Aq ® N 2 )- dimensional max- 
imally entangled states. 

Proof. Since max | (ip \ <f> ) \ 2 is not changed by any 

\4>)e£ 

local unitary transformation on \ip) can always be 
understood in the form of Schmidt decomposition, i.e., 

R-l 

Wab = E I**)- The maximally entangled state can 

i=0 

R-l 

be written as \<j>) = — (U\ <g> U 2 ) \ii), with U\ and 

i— 

U 2 denoting local unitary transformations. Suppose a, = 
(ii |0) , a, G [0, ^], then max \(ip |^)| 2 can be rewritten 



max Utb \d>)\ = max 
\<P)e£ \<P)e£ 
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We arrive at the inequality (7) based on the inequality 

n \ 2 n 

E ) /«< E4 Thus from inequality (8), we can 
i=i / i=i 
obtain 
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(9) 



We would like to emphasize that the '=' in eq. (9) is 
always achieved for the pure states of two qubits. That 
is to say, the right-hand side of eq. (9) is just the con- 
currence of two-qubit pure states. However, for a general 
high dimensional quantum systems, the '=' holds only 
for maximally entangled states (in this case '=' in eq. 
(7) holds), which shows that the inequality (9) provides 
a lower bound for concurrence in high dimension. □ 

Since the maximum of eq. (6) is obtained with \a\ = 

R-l 



1 

>/5' 



can be conveniently chosen as 



i=0 



Therefore, a lower bound of concurrence can be given by 



w>)>v^HH*)(*l)4) (10) 



with p = \ip) (ip\. It is obvious that the lower bound in 
eq. (10) is less than that in eq. (9). But it can be used 
conveniently because there does not exist maximization 
problem. 

The inequality (10) can be immediately gener- 
alized to mixed states. Concurrence for any 
(Aq (g> A^)-dimensional mixed state p is defined as 
C(p) = mmJ2 piC(\(fi)) where the minimum is 
taken over all possible decompositions such that p = 
Eft Ifi) (Vil 7 Eft = 1- Based on the optimal decompo- 
sition p = J2<H \Xi) (Xi\ sucn that C(p) = H c (\Xi)), 
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one can get 
C{ P ) = 

3i 



2R 



R 



2R 
R- 1 



Tr[p 



Tr 



\Xi) (Xi 



(11) 



Inequality (11) holds for any bipartite quantum states 
which provides the key result that will be used later. 



III. EVOLUTION OF CONCURRENCE 

A. One-sided quantum channel 

Next, we will show that eq. (11) can be captured by 
the evolution of some probe states. Let us first consider 
an (N®N) -dimensional bipartite quantum states p with 
only one subsystem undergoing a quantum channel rep- 
resented by the superoperator $i, then the final state 
can be given by pj — where p =Tr[($i ® 1) p] is 

the probability for channel $i which corresponds to non- 
trace-preserving channel [12]. Any quantum state can 
be expanded in a representation spanned by maximally 
entangled states given by 



, N-l 

m - 4= E ^ 



\k)\k®ji),j = Nj +j 1 , (12) 



where jo,ji =0, 1, - ■ ■ ,N — 1, |fc) is the computational 
basis and '©' denotes the addition modulo N. A maxi- 
mally entangled state |$ m ) can also be written as 



|* m > = (QmP- 1 ® 1) \P) 



(13) 



N-l 



where \P) — Y a ij | * J ) , called 'probe quantum state' 

i,j=0 

in this paper, is a generic entangled pure state with full- 
rank P (which can be explicitly written as 



aw 



aoi 
an 



a o(w-i) \ 

Ol(JV-l) 



\ a (N-l)0 «(JV-1)1 ••• a (N-l)(N-l) ) 

or can be directly obtained by the method provided be- 
low eq. (1).) and P -1 denotes the inverse matrix of P. 
$ m are simple unitary transformations determined by eq. 
(12). For example, for the state of a pair of qubits, $ m 
correspond to the three Pauli matrices and the identity, 
respectively. Since any quantum state can be given 

[19], based on maximally entangled state 4>) , by 



VR (V> <g> 1) <f>) = VR (1 ® tp 1 



can also be written as 
4) = (iiP- 1 ® 1) \P) = 



1®1> T (P- 1 ) T )\P), (15) 



where the superscript T denotes transpose. Hence, we 
have 



Tr 



PS 
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= iTr 

V 

1 
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V 



Sp*S{% x ®\)i^){^\, (16) 



where S is the swapping operator defined as S \j) \k) — 

<p) to the 



\j) and we apply eq. (14) and S 

third '='. Eq. (16) can be understood by the Kraus 
representation of superoperator $ x [23]. Based on eq. 
(13) and eq. (15), eq. (16) can arrive at 



Tr 



Pf 



TvSp*S ($i <g> 1) [(l ® 4> T {P- 1 ) 7 ) \P) (P\ 1 ® (P^ 1 )* 



TrSyS m(p- 1 Y 



($10 1) |P) (P| 



(17) 



where p' =Tr($i <g> 1) |P) (P|, p t = p/p' and the star de- 
notes conjugate operation. An alternative derivation can 
be done by substituting p = Y^Pi \<Pi) m t° the first 
line of eq. (16). In a representation of maximally entan- 
gled states, we can have 

p t = Tr [($i <g> 1) p] jp' 

= Ti:^;[($i®l)p]|* m )<$ m |/i/ 



= Tr E 



($i <g> 1) |P) (P| 



Sp*S 



^©(P^f • (18) 



In fact, if the reduced density of the initial state p is 
considered, p t has an alternative and concise form. Let 
Pa =TibPab =TibP with Tr^ denoting trace over sub- 
system B. Consider a decomposition of p = Y^Pi (v»l> 
Pa can be rewritten as = YPififl ■ Thus p t can also 
be given by 



Pt = Tr 



($i <g> 1) |P) (P| 



1® (p-' PA [p-^ 



(19) 

Eq. (19) has a concise form without summation. Sub- 
stitute eq. (18) or eq. (19) into eq. (17), one can find 

that Tr pf <f>^ (ij> has been given by a simple algebra 

(14) on the evolution of the probe state |P) and the original 
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density matrix. That is to say, the lower bound of con- 
currence in cq. (11) can be captured by the evolution of 
the given probe state which can be formally written as 



CI 



2R 
R-l 



Tr 



Bi® l)p] 



R 



, (20) 



where f (x) = 1® {P^f [x] 1® (p- 1 )* S. It 

is obvious that the lower bound of concurrence is deter- 
mined by the evoluted probe state. 



B. Two-sided quantum channel 

Eq. (20) can immediately be generalized to the case 
of two-sided quantum channel, however the form might 
not be as simple as eq. (17). Since the lower bound 
given in eq. (17) is valid for mixed initial states, the 
lower bound for two-sided quantum channel can be easily 
obtained by replacing p in cq. (17) by (1 ® % 2 )p. Thus 

Tr (f)^ ($i ® $ 2 ) p can be written as 

Tr [|<^) ($i ® $ 2 ) = Tr [ ($i ® 1) (1 ® $ 2 ) p 

= ^ Tr E { { m ( p_1 ) T ] ® i) i^> (p\\ [i® (p- 1 ) 



x S[(l®$ 2 )\* pk )(* pk \}*S} 



(21) 



where we have replaced p by a potential decomposition 

of p = l^pfe) (^pk\ which is especially referred to the 

k 

eigenvalue decomposition for simplicity. Applying eq. 
(15) to |* p fe), eq. (21) leads to 



eq. (22) can be rewritten as 

Tr[<A)^ ($i®$ 2 )p 
= ^Tr^($ m |[(l®$ 2 )|P) (P|]|*„). 



rank 



x [(l^p-^pkP- 1 )] 5|$ m )($„|5 
= lTr^($ ro |[(l®$ 2 )|P) (P|]|$„). 



$i®l)|P> (P\T 

(p-h^y 

Bi®l)|P) <P|]* 



(23) 

Finally, it is worth noting that maximally entangled state 
is a special choice of our probe states. Furthermore, the 
value of the lower bound of concurrence does not depend 
on the choice of probe state. It is obvious that eq. (17) 
and eq. (18) correspond to the trace-preserving quantum 
channels. The general results for non- trace-preserving 
channels arc omitted here, which can be directly given 
by adding some normalization constants like the case of 
one-sided quantum channel. 

For integrity, we show that the upper bound given in 
Ref. [8] can be captured by the given probe state |P) , but 
the value of the bound is not changed. Suppose that is 
a bipartite quantum state of qubits, then the concurrence 
can be given by 



C(|V» = 2|det(V)|. 



(24) 



If one of the subsystems undergoes a quantum channel $i, 
the final state can be given by pj = ($i ® 1) \tp) /p\, 
with pi =Tr($i ® 1) \ip) (tp\. Thus we have 



det (V>) 



pip 2 det (P) 



PP ((Jy ® (Ty) p*p (<Ty ® (Ty) , (25) 



Tr^^ ($i®$ 2 )p 

= ^ Tr E { { m ( p_1 ) T ] [($1 ® i) \p) (p\\ [i® (p- 1 )* 

k 

x5(^p- i ®i)*[(i®$ 2 )|p) (p\r ([i , - l ] t *j fc ®i)*s} 
= {[($i ® i) i^> (p\r (i ® p- 1 ^^- 1 ) 

k 

x S [(1 ® $ 2 ) |P) (P|] S (l ® (P-^^P- 1 ) 1 ) } . (22) 



with pp = ($i ® 1) |P) (P| /p2 and 

p 2 =Tr($i ® 1) |P) (P|. Based on eq. (25), we can 
obtain 



C(Pf) = 



det (tp) 



det (P) 



C(PP) " 2|det(P)| ' (26) 



It is obvious for a mixed initial state g that eq. (25) is 
extended to 



C(Qf) < 



C(g) ■ C{ PP ) 



Substituting eq. (21) into eq. (17), one can find that the 
lower bound of concurrence C[($i®$ 2 )p] can be cap- 
tured by the evolutions of the probe state under the two 
quantum channels. In order to avoid the decomposition 
of the initial state p, one can expand [(1 ® $ 2 ) |P) (P|] in 
the representation of maximally entangled states. Thus 



2|det(P)| ■ 

For two-sided quantum channel, one can easily obtain 
C(g)-C(p P1 ) C(p P2 ) 



(27) 



C{Qf) 



< 



2|det(P)| 2|det(P)|' 



(28) 



with PP1 ($! ® 1) |P) (P| /p 2 , PP2 

(1 ® $ 2 ) |P) (P| /p' 2 and p' 2 =Tr(l ® $ 2 ) |P) (P| . 
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FIG. 1: (Dimensionless)The upper bound of concurrence 
(dashed line) given in Ref. [12], the concurrence (solid line) 
and the lower bound (dash-dot line) of quantum state pf vs 
x. pf is the final output state of an initial state p with the 
subsystems undergoing a quantum channel $i and $2, respec- 
tively. 



Thus concurrence (especially bipartite concurrence of 
qubits) can be better evaluated by the lower bound and 
the upper bound given in eq. (28) than by only one 
bound. From cq. (20) as well as eq. (22), one might 
think that it is not so convenient compared with the up- 
per bound in Ref. [8] , because the upper bound is a sim- 
ple linear relationship between the concurrence of evolved 
probe state and that of the initial states. However, gen- 
erally speaking, concurrence per se is not a direct ob- 
servable [24-26], therefore, to evaluate concurrence of a 
state in practice, one has to evaluate the quantum state 
by quantum state tomography [27] and then turns to a 
mathematical procedure. In other words, it is inevitable 
for Ref. [8] to evaluate quantum states in practical sce- 
nario. In this sense, we think that their practicability is 
almost the same. What is more, our lower bound has 
obvious advantages: 1) The lower bound has the consis- 
tent spirit with entanglement measures of mixed states 
for which the lower bounds (infimum) are usually needed; 

2) Bounds, especially lower bounds, with elegant forms 
like Ref. [8] might be difficult to provided. In particular, 
so far there have not been analytic results of entangle- 
ment measure (in particular, concurrence included) for 
general high-dimensional quantum systems, therefore, if 
the bounds of entanglement for high-dimensional mixed 
states still include the calculation of high dimensional 
entanglement measures, it only formally provides an ele- 
gant relationship, but it has usually poor practicability. 

3) The derivation based on our extended lower bound 
provides a universal method for all analytic bounds of 
entanglement measure, with which one can only focus all 
the attention on the tightness of the bounds, but there 
might be a great deal of difference on the complexity of 



the final results between different lower bounds. 



IV. A SIMPLE APPLICATION 

As an application, we only consider a bipartite sys- 
tem of qubits, because there exist analytic concurrence 
and acceptable upper bounds for bipartite mixed state 
of qubits. Thus one can directly find the tightness of 
our lower bound by comparing it with the concurrence 
and the upper bound. The bipartite quantum state we 
considered is given by 

(1 - x) 

p = x Pr + K - ' 1, x €[0,1] (29) 

with 

/ 0.4322 0.2113 0.1073 0.3369 \ 
_ I 0.2113 0.1845 0.0406 0.1798 I , > 

Pr ~ I 0.1073 0.0406 0.0504 0.1144 ' 
V 0.3369 0.1798 0.1144 0.3330 / 

randomly generated by Matlab 6.5. We suppose that 
each subsystem of p undergoes an amplitude-damping 
quantum channel given in Kraus representation [28] as 




The final state can be written as 

p/ = ($i®$ 2 )p. (33) 

The upper bound, the concurrence itself and the lower 
bound are given in FIG. 1, respectively, from which we 
can find that our lower bound is a good evaluation of 
concurrence. The probe state can be chosen freely. 



V. CONCLUSION AND DISCUSSION 

We have presented a lower bound of concurrence. In 
particular, when the quantum state under consideration 
evolves under a quantum channel, the lower bound of 
concurrence can be completely captured by the evolution 
of probe states. Thus, the evolution of concurrence for 
any initial state can be well evaluated by the lower bound 
and the upper bound. Furthermore, the lower bound is 
also suitable for high-dimensional quantum state. We 
would like to emphasize that, even though the form of 
the lower bound seems not to be as elegant as that of the 
upper bound, but their practicality is almost the same. 
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